Astronomy & Astrophysics manuscript no. ForarXiv 


February 9, 2012 


(DOI: will be inserted by hand later) 





The orbit of 2010 TK7. Possible regions of stability for other Earth 

Trojan asteroids 

R. Dvorak^, C. Lhotka^, L. Zhou^ 



O 
<N 

00 



o 

03 



(N 
> 

00 

o 

<N 

o 

(N 



^ Universitats Stern warte Wien, Tiirkenschanzstr. 17, A- 11 80 Wien, Austria, 

^ Department de Mathematique (naXys), Rempart de la Vierge, 8, B-5000 Namur, Belgium, 

^ Department of Astronomy & Key Laboratory of Modern Astronomy and Astrophysics in Ministry of Education, Nanjing 
University, Nanjing 210093, China 

Received; accepted 

Abstract. Recently the first Earth Trojan has been observed (Mainzer et al., ApJ 731) and found to be on an interesting orbit 
close to the Lagrange point L4 (Connors et al.. Nature 475). In the present study we therefore perform a detailed investigation 
on the stability of its orbit and moreover extend the study to give an idea of the probability to find additional Earth-Trojans. Our 
results are derived using diff'erent approaches: a) we derive an analytical mapping in the spatial elliptic restricted three-body 
problem to find the phase space structure of the dynamical problem. We explore the stability of the asteroid in the context of the 
phase space geometry, including the indirect influence of the additional planets of our Solar system, b) We use precise numerical 
methods to integrate the orbit forward and backward in time in diff'erent dynamical models. Based on a set of 400 clone orbits 
we derive the probability of capture and escape of the Earth Trojan asteroids 2010 TK7. c) To this end we perform an extensive 
numerical investigation of the stability region of the Earth's Lagrangian points. We present a detailed parameter study in the 
regime of possible stable tadpole and horseshoe orbits of additional Earth-Trojans, i.e. with respect to the semi-major axes and 
inclinations of thousands of fictitious Trojans. All three approaches underline that the Earth Trojan asteroid 2010 TK7 finds 
himself in an unstable region on the edge of a stable zone; additional Earth-Trojan asteroids may be found in this regime of 
stability. 



1. Introduction 

The giant planets Jupiter and Neptune are known to host Trojan 
asteroids, and also Mars (Bowell et al. 1990) is hosting several 
co-orbiting asteroids. These kind of bodies move in the same 
orbit as the planets, but around 60° ahead or 60° behind the 
planet close to the so-called Lagrange points L4 or L5 (see 
Fig. [B. It was a great surprise when q uite recently the firs t 
Earth Trojan 2010 TK7 was discovered ("Mainzer et al. 201 1"). 
Although many studies have shown the principal possibility of 
their existence until this event all attempts of finding one were 
unsuccessful. Another small asteroid in the 1:1 mean motion 
resonance (MMR) with the Earth was found earlier in 1986, by 
Duncan Waldron and this asteroid (3753 Cruithne) was later 
identified as a celestial body in a horseshoe orbit around both 
equilateral Lagrange points of the Earth being thus not a 'real' 
Trojan in the sense of its original definition (the same is true 
for the recently found asteroid 2010 SO 15) . 

Many theoretical studies exist to establish the stability 
of t he Lagrange points in simplified ni o dels, i.e. the works 
of dRabe 19671: leien & Schubart 1984 iLhotka et al. 2008 : 



Erdi et al. 2QQ9h , beyond many others. Ever since extensive nu- 
merical studies have been undertaken in finding the extension 
of the stability regions around the equilibrium points for the 
planets, e.g. in (Schwarz et al. 2004; "Dvorak & Schwarz 2005|; 
Robutel etal. 20051: ^reistetter 2006; Dvorak et al. 2007J. 
Especially importa nt results are due to the work of 
dMikkola & Innanen 1992: iTabachnik & Evans 20o3; 

Brasser & L etho 2 00^ a nd dSchol l & Marzari 2004). In 
( Mikk ola & Innanenl992h and dZhan g & Innanen 1995) the 
authors found that Venus, Earth and Mars can host co-orbital 
asteroids up to 10 Myr s. According to the investigation of 
dTabachnik & Evans 2000) Earth's Trojans are on stable orbits 
when their inclinations are relatively low (/ < 16°); a second 
stability window exists according to them for 16° < / < 24°. 

Morals & Morbidelli (2002) studied the orbital distribution 
of the near-Earth asteroids (NEAs) that experience the episode 
of being in the 1 : 1 MMR w ith the Earth. In a most recent study 
dSchwarz & Dvorak 201 ih the possibility of captures of aster- 
oids by the terrestrial planets into the 1:1 MMR was investi- 
gated and many temp orary captures including jumping Trojans 
dTsiganis et al. 2000b were found; as we will see also 2010 TK7 
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has such an interesting captured orbit. 

In our investigations - initiated by the finding of 2010 
TK7 - we concentrate on Earth Trojans. Three diflferent ap- 
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proaches should clarify the stability problem of these aster- 
oids: in Chapter 2 we perform analytical studies in a simpli- 
fied model, valid on short time scales. We investigate the phase 
phase structure and the influence of additional perturbations in 
the simplified model. In Chapter 3 we study in great detail the 
actual orbit of 2010 TK7 together with 400 clone orbits to ob- 
tain a better statistics to state the probability of capture and es- 
cape of this Trojan asteroid. To this end we perform in Chapter 
4 an extensive numerical investigation of the parameter space 
relevant to determine the stable and unstable regions close to 
the equilibrium point L4 and show the interplay between sec- 
ular resonances and the stability of motion. The summary and 
conclusions of the present study are found in Section 5. 




Fig. 1. Geometry of the circular restricted three-body problem 
in a rotating coordinate system: mass of the Sun jiq and the 
Earth //'; the unstable equilibria ^1,2,3 are located on the x-axis; 
the stable points L4 and L5 are found 60° ahead and behind 
the primary mass jd\ The color code defines the areas of the 
equipotential. In the elliptic problem the rotating coordinate 
system has to replaced by a non uniformly rotating and pul- 
sating reference frame. 



2. A symplectic mapping model 

The Hamiltonian describing the motion of 2010 TK7 in the spa- 
tial elliptic restricted three-body problem (Sun-Earth-asteroid) 
takes the form: 

H = i^Kep + 7^ + fi'R(a, e, U CO, Q, M, M'\ p') (1) 

where Hxep denotes the Keplerian part, the function [I'R is the 
perturbation due to the Earth with mass jd' and the variable T 
is the action conjugated to time (assuming, that the mean mo- 
tion n' of the Earth is equal to one). Moreover, a, e, /, co, Q, M 
are the semi-major axis, the eccentricity, the inclination, the 
perihelion, the longitude of the ascending node and the mean 

^ both Lagrange points have in principal the same dy- 
namical beh aviour a s we know from ea rlier studies e.g. 
(iNesvornv &Vokrouhlickv 20091 : IZhou et al. 2009i) 




60 120 180 240 300 

T=A-A' 



Fig. 2. Phase portrait - projected to the (r, r)-plane. The Earth 
is located at (0, 1), the fixed points L4, L3, L5 (corresponding to 
the equilibria of the averaged system) are situated at (60°, 1), 
(180°, 1) and (300°, 1), respectively. The projection of the mean 
orbit of the asteroid 2010 TK7 for ±6000 years is shown in red. 

anomaly of the asteroid, respectively, while M' denotes the 
mean anomaly of the Earth. 

In contrast to the classical expansion techniques of the per- 
turbing function R we do not replace the remaining orbital pa- 
rameters of the Earth p' = (a\ e\ i\ co' ,Q.') with their numer- 
ical values but rather keep them as free parameters in the on- 
going calculations. The simple reason for that is to be able to 
investigate their influence on the dynamics of the massless as- 
teroid with time, i.e. to see their influence on the phase space 
geometry, later on. For the mapping approach we set the sys- 
tem of constants (gravitational constant G and total mass of 
the system) equal unity such that G(ji{) +11') = 1 which also 
implies that a' as well as a are close to unity and in these 
dimensionless units one revolution period of the Earth takes 
P = In. Furthermore, we use as action-angle variables the 
modified Delaunay variables which are defined in terms of the 
classical Delaunay variables by: 

/ = M, g = 61;, /z = Q 

L= ^fa, G = Lil-e^, H = Gcos(0 

as 

Al = I -\- g -\- h, A2 = -g - h, A3 = -h 
Al = L, A2 = L - G, A3 = G - if 

(and similar for the primed variables of the Earth). In this set- 
ting, the Hamiltonian ([T]) becomes: 

H=--^ + T+fi'RiA,KA[;q') (2) 

ZA2 

where we used the short-hand notation A = (/li, /I2, /I3) and 
A = (Al, A2, A3). In addition we abbreviate the parameter vec- 
tor q' = (/l2,/l3, Ap A2, to define the vector of modified 



Dvorak et al.: The Earth Trojan asteroid 2010 TK7 



3 



Delaunay variables for the Earth. Next, we implement a sym- 
plectic change of coordinates and momenta suitable to describe 
the motion of the asteroid close to the 1 : 1 MMR. We define 
the resonant angle and conjugated momenta as: 

T = Ai-A[, T = Ai (3) 

while the other variables transform by the identity: 

M = (f, A3 = 6, A2 = O, A3 = 

and T' = Al - A[. We replace Ai in Q according to ^ and 
implement the standard averaging procedure over the fast angle 
A' via the formula: 



H = -^^^jy'R (A, A, A[ ; q') dA' , 



(4) 



with P = 271, to get the averaged Hamiltonian function of the 
form H = H(i//, ^; q') with i// = (r, cp, 6) and ^ = (7, 0, 0). Our 
aim is to construct a symplectic mapping which transforms the 
state vector = (jk^ Ok, Tk, Oy^, &k) at discrete times k 

(multiples of P) to the state vector (i//k+i,^k+i) at times ^ + 1. 
For this reason we define the generating function 



q' - Wq' in, (fk 



q' y^k^ ^k^ ^k, Tk+l, ^k+l, 0)^+1 ; ^0 - 

i//k ' Jk+i+ ^7iH {if/k, Jk+i ; q') 



where the symbo l • is the dot pr oduct. A s it h as been s hown and 
already used in (|Had iidemetri ou 1992 ;^ Had iidemetriou 19991: 
Hadiidemetriou & Vovatzis 20 00) or (Ferraz-M ello 1997h the 
generating function defines a mapping of the form: 



Tk+l 



dTi 
dTk 



k+i 



(fk+i 



0,= 



^k+i 



'k+i 



difk 



&k = 



(5) 



System © defines a symplectic change of coordinates on 
the Poincare surface of section obtained by the averaged sys- 
tem defined by ©. The set of variables i//) is related to the 
mean orbital elements of the asteroid 2010 TK7 via another 
generating function (not derived here) which defines the aver- 
aging process given in (|4|. The system © therefore describes 
the evolution of the mean orbital elements of the asteroid at 
discrete times t = k- P. Note, that the system is implicit to pre- 
serve the Hamiltonian structure of the original problem. For 
given initial conditions (i//k^o,^k^o) it can either be iterated 
by solving the system of diff'erence equations implicitly for 
C¥k+i, ^k+i) or by ni ore sophisticated procedures as described 
e.g. in dLhotka 20091) . 

A typical projection of the phase portrait to the (r, r)-plane 
is shown in Fig.O The lines in black were obtained by varying 
T within (0°, 360°) along T = I. The plot also shows the mean 
orbit of 2010 TK7 in red. It has to be compared with Figure 2 of 
dConnors et al. 201lk : while the authors derive the red curve by 
numerical averaging, the averaged orbit of the present approach 
is based on Eq. (O. 

As it is well-known, e.g. (iMorais & Morbidelli 2OO2I: 
Brasser & Letho 20021: ISchoU et al. 2005*). the influence of the 
other Solar system bodies, i.e. the direct influence of the major 



planets, afl'ects the motion of the asteroids on long time scales. 
The mapping does not take into account these direct efl'ects. 
Thus, the mean orbit as shown in Fig. [2l can only been seen 
as a first approximative solution of the problem. Moreover, the 
simulations were done using low order expansions of the per- 
turbing function, since we wanted to keep the dependency of 
the model on the orbital parameters of the Earth. The resulting 
mapping model is therefore only valid within a good conver- 
gence regime of the Fourier-Taylor series expansions used to 
approximate 7? in ([T]). In addition, since the orbital parameters 
of the asteroid 2010 TK7 may reach high values (as found from 
numerical simulations) the error of the approximation of the 
series expansions may exceed the cumulative eff'ect of the per- 
turbation eff'ects due to the Solar system and may also change 
the picture dramatically on longer time scales. An estimate of 
the error terms and higher order series expansions together with 
their influence on the long term dynamics of the asteroid 2010 
TK7 is currently under investigation. 

3. The orbit of 2010 TK7 

The Earth Trojan 2010 TK7 fundamentally is a NEA, and its 
orbital elements (listed in Table [B can be found also on the 
AstDyS (Asteroids - Dynamic Site) website 0. At first glance, 
this object has a large eccentricity (~ 0.19) and it is reasonable 
to suspect that it's on a unstable orbit. To check the orbital sta- 
bility of this object, we perform numerical simulations of the 
orbit. 

We adopt two dynamical models in our simulations. Both 
models contain the Sun and 8 planets from Mercury to 
Neptune. They diff'er from each other by diff'erent settings of 
the Earth-Moon system. In one model, the Earth-Moon sys- 
tem is simply treated as one mass point with the combined 
mass of the Earth and the Moon at the barycenter of the sys- 
tem. And in the other model, the Earth and the Moon are re- 
garded as two separated objects as in reality. Hereafter we de- 
note the former by EMB model and the latter by E-hM model 
for short. By comparing these two models, we may testify the 
reliability of the mod el adopted later in this paper and also by 
(IConnors et al. 201V ) in their numerical simulations. 

Considering the uncertainties in observation and orbital de- 
termination, it is necessary to study clone orbits within the error 
bars. Beside the nominal orbit, 400 clone orbits are generated 
using the covariance matrix given by the AstDyS website. As 
listed in Tabled the errors are very small, therefore the initial 
conditions of the clone orbits are very close to the nominal orbit 
as shown in Fig.[3l We show here the distribution of a, e, i and 
II, the other two elements M, co not illustrated here have sim- 
ilar distributions according to the corresponding uncertainties 
listed in Table [T] 

Starting from the initial conditions generated in the 
above way, we integrate their orbits up to 1 million years 
(Myr) in both directions (forward to the future and back- 
ward to the past ). We used the integrator package Mercury 6 
dChambers I999I) and verified some results using the Lie- 
integrator (IHanslmeier & Dvorak 19841) . The comparison be- 
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Table 1. The Keplerian elements of 2010 TK7 given at epoch 
JD2455 800.5. The elements and I - cr variations are from the 
AstDyS (see text). The covariance matrix is also taken from the 
same website. 





Value 


1-cr variation 


Unit 


a 


1.00037 


2.546 X 10"/ 


AU 


e 


0.190818 


9.057 X 10-^ 




i 


20.88 


7.274 X 10-5 


deg. 


Q 


96.539 


1.842 X 10-4 


deg. 


CO 


45.846 


2.309 X 10-4 


deg. 


M 


217.329 


1.848 X 10-4 


deg. 



tween the results from different integrators shows that the re- 
sults are consistent with each other. 



20.878 20.879 20.880 20.881 1.000355 1.000360 1.000365 1.000370 1.000375 

i (degree) a (AU) 

Fig. 3. The initial conditions of clone orbits. Dots are for 400 
clones, and the nominal orbit is indicated by the red star. We 
show here in two panels only the inclination, ascending node, 
semimajor axis and eccentricity. 

For an asteroid on a Trojan-like orbit, the critical angle (res- 
onant angle) is the difference between the mean longitudes of 
the asteroid A and the corresponding host planet (here it's the 
Earth-Moon system) X\ like the r in Eq. ([3]). In the EMB model 
the resonant angle is /I - /Iemb , while in the E-hM model we cal- 
culate the position and velocity of the barycenter of the Earth- 
Moon system from the orbital elements of the Earth and the 
Moon, and then compute the resonant angle in the same way as 
in the EMB model. When the resonant angle A- X = A- Aemb 
librates around 60°, i.e. 0° < A- A' < 180°, the asteroid is said 
to be an L4 Trojan, when it librates around -60° or 300°, i.e. 
180° < /I - /I' < 360° it is an L5 Trojan, when the asteroid 
librates with an amplitude larger than 1 80° it is on a horseshoe 
orbit, and finally when A - A' circulates the asteroid leaves the 
Trojan-like orbit. 

In Fig.m we illustrate the temporal evolution of the nomi- 
nal orbit in the E-hM model. The first impression obtained from 
its behaviour may be that it's a temporal Earth Trojan, i.e. judg- 
ing from the resonant angle it was not an Earth Trojan 0.055 
Myr before and it will not stay on a Trojan-like orbit after 0.37 
Myr. The period of its being an L4 Trojan is even much shorter, 
less than 2,000 years in the past and less than 17,000 years 
in the future (as partly shown in Fig.O. During the period of 
being a Trojan, the semimajor axis a shows regular variations 
librating around 1.0 AU with a small amplitude, however, the 



variations of eccentricity e and inclination / reveal the chaotic 
character of the motion, as shown in Fig.|4l 




-0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6-0.2 -0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 



Time (Myr) Time (Myr) 

Fig. 4. The temporal evolution of the nominal orbit in 1 million 
years, both forward and backward. In four panels we show the 
evolutions of the semimajor axis (upper left), eccentricity (up- 
per right), inclination (lower left) and the resonant angle (lower 
right). 

The chaotic character of the motion implies not only a sen- 
sitive dependence on the initial conditions but also a sensi- 
tive dependence on the dynamical model. To compare the mo- 
tions in two models, i.e. the EMB model and E-hM model, we 
show as an example in Fig. [5] the evolution of the resonant an- 
gle of the nominal orbit in both models. Around the starting 
point, two curves representing the motions in two models are 
almost the same so that they overlap each other exactly. But 
the difference between them becomes distinguishable only af- 
ter about 2,000 years in both directions. This difference on one 
hand arises from the different settings of the models, on the 
other hand, it is due to the chaotic character of the motion. In 
this sense, it is impossible to draw any convincible conclusion 
about these two dynamical models on long time scales by com- 
paring just the motions of an individual orbit in two models. 
So we turn to analyze statistically the four hundred clone orbits 
below. 

As we mentioned above, 400 clone orbits are calculated in 
both models and in both directions of time. At starting moment 
^ = these clones are around the L4 Lagrange point as the as- 
teroid 2010 TK7, but in the evolution, the objects may jump 
from the L4 region to L5 region, or they may move from tad- 
pole orbit to horseshoe orbit, and even, they may escape from 
the Trojan-like orbit (the 1:1 MMR). To examine the motion, 
we check the resonant angle A - Aemb at each step during our 
simulations. When A - Aemb for the first time is larger than 
180°, the object is regarded as leaving the L4 region, and we 
denote the time at this moment as ti . When A - Aemb reaches 
for the first time 360°, the object escapes the 1 : 1 MMR, and we 
denote this moment as t2. 

The ti and t2 of 400 clone orbits are summarized in Fig. [6] 
and Fig. [71 For ti, most of the clones will leave the L4 region in 
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-5000 -4000 -3000 -2000 -1000 1000 2000 3000 4000 5000 

Time (yr) 

Fig. 5. The resonant angle of the nominal orbit in two models. 
The difference between them is distinguishable very soon (~ 
2000 years) mainly due to the chaotic character of the motion 
(see discussion in text). 



- •□- • Backward EMB 
—•—Forward EMB 
V77A Backward E+M 
■■■d Forward E+M 




log(^i) 



Fig. 6. The distribution of the time when clone orbits leave 
from current L4 region. The escape time from the tadpole re- 
gion is given in years in logarithm scale. 

about 1.5 X 10"^ years in future, while the backward integrations 
indicate that the clones "entered" the L4 region in about 1.8 x 
10^ years ago (Fig.O. 

In the backward integrations, the earliest escape from the L4 
region happens at ^ ~ 1440 years both in the E+M and EMB 
models. But there are only a few such orbits. Most of them es- 
cape in the time range between t ~ 1680 and t ~ 1820 years 
which causes the peak centered at log(^i) = 3.25 in Fig.O In 
fact, from Fig.[5l we may derive the libration period of the reso- 
nant angle is around 350 year^, and we also note that the shift 
of the orbit from L4 to L5 region happens when the resonant 
angle reaches the maximum in a libration period. The libration 
period and amplitude must be tuned by other periodic effects 
(e.g. secular resonances). After a number of complete periods 
of evolution, the resonant angles of some clone orbits reach 



^ Due to the large libration amplitude, this value is nearly 50 per- 
cent larger than the synodic period of a tadpole orbit which can b e es- 
timated using the formula (27yL//4)"^/2 ( Murray & Dermott 199^ . For 
an Earth Trojan, the mass ratio \i - 3.04 x 10~^ and the period is 
~ 220 years. 



their maxima, and they do not librate back but escape from 
there towards the neighbourhood of L5 . After another complete 
libration period, many more clone orbits escape in the same 
way as before. That's why the escapes from the L4 region seem 
to happen more or less suddenly (at t ~ 1440 -h 350 = 1790 
^q3.25 years). All clone orbits are in a tiny region confined by 
the error bars, so that they suffer nearly the same dynamical ef- 
fects in a short timespan. Eor those clones staying in the region 
longer, they spread in the phase space and thus suffer different 
dynamical effects and consequently their escape times diverge. 

The distributions of t\ for the EMB model and the E-hM 
model in Fig. [6] match each other very well, therefore we may 
draw another conclusion here that the difference between these 
two dynamical models is ignorable, or in other words, the sim- 
plified EMB model is a reasonable and reliable model for in- 
vestigating the Earth Trojans' long term dynamics. Particularly, 
we would like to note that the asteroid 2010 TK7 and the clones 
all have very large libration amplitudes around the L4 and/or 
L5 point (see for example Fig.O, i.e. they may approach to the 
Earth-Moon system (mean longitude difference between them 
may be less than 20°). But in fact, even when an Earth Trojan 
approaches the Earth to within 10° apart in mean longitude, the 
distance between the asteroid and the Earth is about 0.18 AU, 
which is still about 70 times larger than the distance between 
the Earth and the Moon. Taking into account the inclination of 
the Earth Trojan's orbit, the distance between the Trojan and 
the Earth is even larger. 



- □- Backward EMB 
— •— Forward EMB 
V77A Backward E+M 
Forward E+M 




iog(g 

Fig. 7. The distribution of the time when clone orbits escape 
from the 1 : 1 MMR, both in forward and backward integrations. 
The escape time is given in years in logarithm scale. 

As for t2 distribution in Fig. [71 we see the peak of escaping 
time in the backward integrations locates at ~ 3.0 x 10"^ years. 
As we mentioned above, the Earth Trojan 2010 TK7 should be a 
temporal Trojan. The capture of this asteroid onto the 1 : 1 MMR 
with the Earth happened most probably 30,000 years ago. In 
the forward integrations, most of clones will escape from the 
MMR around 2.0 x 10^ years after. The total lifetime of this 
asteroid in the 1:1 MMR is less than 0.25 Myr. And again, the 
two models E-hM and EMB give nearly the same distributions. 

In our calculations, the episode of clones being L4 Earth 
Trojans lasts typically a little less than - 17, 000 years. As for 
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the lifetime of clones staying in the 1:1 MMR, there are 10 
clones surviving 1 Myr in the backward integrations in both 
models, while for the forward integrations, 12 in the EMB 
model and 19 in the E+M model stay in the MMR till the end of 
integrations. But none of the clones survives in both temporal 
directions. In one word, less than 5 percent of clones stay in the 
MMR up to 1 Myr. Morals & Morbidelli (2002) calculated the 
probability of a NEA being captured into the 1 : 1 MMR ("coor- 
bital orbit" in their paper), and they found that each episode of a 
NEA being coorbital on average is 25,000 years and none lasts 
longer than 1 Myr. Our result does not conflict with their con- 
clusion, because our calculations are for the individual asteroid 
2010 TK7, and its eccentricity (~ 0.2) is smaller than the typi- 
cal eccentricity in their samples (most of them have e > 0.28). 

Some fluctuations in the distribution of ti and t2 (by "fluc- 
tuation" we mean more peaks deviated from a trivial normal 
distribution) can be found in both models in Fig. Sand Fig. [71 
In these figures, the time is given in logarithm to include a wide 
time range. If we plot the time linearly, these fluctuations show 
some periodic character, implying that some periodic mecha- 
nisms (e.g. secular resonances) are afl'ecting the libration am- 
plitudes. However to depict a secular resonance map is a long 
story, we would like leave it to a separate paper but show first 
results in Section 4.3. 



4. Determination of the stable regions 

As was pointed out in the previous sections the orbit of 2010 
TK7 is not very stable. The goal of this part of our investiga- 
tion is to determine the largeness of the stable regions arou 
the stable equilibrium points. This can be done in a real 
tic dynamical model only with the aid of numerical integi 
tions where many fictitious Trojans are checked concerni 
the stability of their orbits. The dynamical model used wae 
'truncated' planetary system with the planets Venus to Satu 
(Ve2Sa). The simple reason for not taking into account all t 
planets is that it would at least require four times longer CI 
times on our computers. It is mainly due to the orbit of the i 
nermost planet Mercury which demands for a step size of o 
quarter of the one of Venufl Test computation for selected si 
ble and unstable orbits showed that the qualitative behave 
of such orbits is not a diff'erent one and consequently also t 
largeness of the stable region can be regarded as the 'real' 01 
It means that the gravitational perturbation of Mercury on 
Trojan for this study is neglectable. The Moon was not e 
plicitly integrated but the barycenter Earth Moon with the i 
spective masses was taken as one body. As has been shown 
Chapter 3 there are only very small diff'erences when one i 
tegrates the Trojan orbits in these models, on long time scales. 
On the contrary Jupiter's direct perturbation on an Earth Trojan 
is very large and cannot be neglected especially for large libra- 
tions which bring the Trojan close to the Lagrange point L3. 

The integration method was the Lie- series method which 
we used already quite often in similar inves tigations (e.g. for 
Neptune ( Zhou et al. 20091: Izhou et al. 201 ih and recently also 



Uranus' Trojans ('Dvorak et al. 20 id) ). This method is b ased 
on a work by (Hanslmeier & Dvorak 1984 iDelva 19841) and 
(ILichtenegger 1 984) and was slightly modified for our problem 
(see also ("Eggl & Dvorak 201 0)). It has an automatic step-size 
control and turned out to be fast and precise. 

The initial conditions were chosen in the following way: 
for the fictitious Trojans the orbital elements M, Q, e were set 
to the one of the Earth. The perihelion was always oJEarth + 60° 
0, the semimajor axes for the Trojan was set to slightly smaller 
and larger values to cover the stable region (along the connect- 
ing line between Sun and the Lagrange point). We started with 
integrations in the plane of the orbit of the Earth and changed 
in additional runs the inclinations up to / = 60°. 

As to distinguish between stable and unstable orbits we 
used diff'erent indicators: the most straightforward was to check 
the eccentricity of the object because there turned out to be 
a sharp cut for ^ > 0.3: any orbit which achieved this value 
during the integration left the area around the Lagrange point 
(checked by its distance to this point). In another test we com- 
puted directly the escape times and finally we also computed 
the libration width of the Trojan which is a well established 
check of stability for a Trojan. The length of integration was 
thoroughly chosen; although we had a computer grid available 
the main computations covered only 10^ years but some tests 
has been undertaken up to 10^ years. 

4.1. The stable regions for different time scales 




semimajor axes 



^ The inclusion of Uranus and Neptune would not dramatically 
change this integration time 



Fig. 8. Stability diagram using emax for the window between 
/ < 60 (y-axes) for diff'erent initial semimajor axes (x-axes) of 
the fictitious Trojans; the color stands for the maximum eccen- 
tricity during the integration of 10^ years. 

The stability diagram for a cut through the Lagrange point 
L4 is shown in Fig. [8] where we plotted the maximum eccen- 



which means that the position was the Lagrange point L4 
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tricity of the orbits of 100 massless fictitious Trojans for d 
ferent values of the semimajor axes (x-axes) versus the inc 
nations of the Trojan (y-axes). One can see that for small i 
clinations the stable region extends in semimajor axes betwe 
0.99 < a < 1.01, then there is a decrease of the largeness 
the region visible up to about / = 20°. Inside this region i\ 
finger-like slightly less stable regions are visible on both sid 
of the a V-shaped structure centered at a = I leading frc 
15° < / < 20°. Then an unstable strip follows to about / = 2 
with some connection to a larger rectangular stable region wi 
0.997 < a < 1.003 and 26° < / < 40°. Stable symmet 
fingers inside this region are clearly visible and build a conti 
uation of the unstable fingers around a ^ 0.997 and a ^ 1.0 
for the large stable area for / < 20°. For larger inclinatio 
on / > 40° the unstable regions extends for all values of t 
semimajor axes with the exception of a small U-shaped stal 
window around / = 50°. Note that inside this unstable regi 
(red to yellow) sometimes small stable island seem to appe 
on both sides of the stable region / < 40° which - after a lon^ 
integration - disappear. 

Concerning the structures visible in the figure men- 
tioned above we note that similar ones have been dis- 
covered for the Trojan regions of the outer p lanets 
in the papers by (IMichtchenko & Ferraz-Mello 2002h and 



(iNes vorny & Po nes 2002b for larger eccentricities. In our pa- 
per we investigated the stability for larger inclinations; the dif- 
ferent features inside stable r egions for Earth T rojans are com- 
parable to the ones found by fehou et al. 2009") for the Neptune 
Trojan region. The instabilities inside the stable region are 
caused by secular resonan ces - as it has been show n in detail 
for the Jupiter Trojans by (iRobutel & Gabern 2006"). In our ar- 
ticle we just want to show - for the moment - the complicated 
structures. To make a detai led analysis conc erning the reso- 
nances like in the paper by ('Zhou et al. 2009*) or even the for- 
mer mentioned article by ( Robutel & Gabern 2006) takes quite 
a long time; we started this work already. 

In a next step we extended the integration to 10^ years in 
the same dynamical model Ve2Sa. The region for small incli- 
nations is still stable with the same extension in semimajor axes 
(Fig.O. New features, unstable vertical strips, appear for dif- 
ferent values of a and increasing values of the inclination. For 
a ^ 0.997 and a ^ 1.003 these regions were already visible as 
being less stable in the former plot (Fig. [8]) for small inclina- 
tions. New unstable vertical fingers appear also for a ~ 0.995 
and a ~ 1.005 and 10° < / < 16°. A new characteristic is the 
appearance of unstable fingers for a ~ 0.998 and a ~ 1.002 
with / < 6°. The V-shaped unstable region centered at a = 1 
was already visible in Fig.[8]but the detailed structures become 
more clear using a longer integration time. 

In the next Fig[TOl the stable window for 28° < / < 40° 
was studied separately (also for 10^ years). Again one can see 
that the most stable regions visible through the dark blue - in- 
dicating small e-max values - are close to the edges in a. These 
edges are quite irregular and some smaller unstable islands ap- 
pear inside; the whole stable region is somewhat tattered. What 
we do not show here is the disappearance of the U-shaped is- 
land for / ~ 50° for a longer integration time. 



12 
10 



0.998 1 1.002 1.004 1 .C 
semimajor axes 



Fig. 9. Caption like for ¥ig.\M 



We note th at we disagree, in some aspects, with the re- 
sults given by (Tabachnik & Evans 20 00) who claim that sta- 
ble regions for Earth Trojans are possible for 24° < / < 
34°; the stable window we found is shifted outwards to 
larger inclinations. We also disagree with the later paper by 
(IBrasser & Letho 20021) where they say that orbits with an in- 
clination of 12° < i < 25° on are unstable. In the last work the 
authors quite nicely determine the secular frequencies involved 
leading to unstable motion. We already started additional com- 
putations where the determination of the resonances are un- 
dertaken, but their detailed analysis will be shown in a longer 
article (in preparation). 



4.2. The libration amplitudes 

To determine the libration width of the region we choose for 
the semimajor axes a grid of 0.995 < a < 1.005 for 50 diff'erent 
fictitious bodies equally distributed in the mentioned interval 
where the inclination was set to values 0° < / < 56°. The inte- 
gration time was only 10^ years for this study. The colors (from 
blue to yellow) indicate the amplitude of libration; we can see 
a well defined stable region in the range 0.997 < a < 1.003 
for inclinations i < 19°. This rectangular like region seems to 
contradict the results shown in Fig. [HI but a closer look shows 
- through the color of dark yellow - that the libration angle is 
in the order of 170°; thus on the edge we have horseshoe orbits 
which enclose both equilateral Lagrange points. An unstable 
strip then arises 20° < / < 24° where close to a - 1.0005 and 
/ ~ 21° the 2010 TK7 is located. Then we see again a large 
almost rectangle like stable region 25° < / < 39° for the same 
range of a as the first stable region of tadpole orbits for only 
slightly inclined Trojan orbits. Note that the edges are sharp 
and no horseshoe orbits are located there. From / = 40° no 
more stable orbits exist with the exception of a small stable 
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Fig. 10. Stability diagram using emax for the window between 
26 < i < 40 (y-axes) for different initial semimajor axes of the 
fictitious Trojan; the color stands for the maximum eccentricity 
during the integration of 10^ years. 



window with tadpole orbits for / ~ 50°; these orbits will be 
unstable for longer integration time as indicated above. 



tially being in a horseshoe orbit (Fig. [121 lower graphs). On the 
left pictures of this figure the semimajor axes show a periodic 
change between a maximum and a minimum value close to the 
semimajor axis of the Earth for the tadpole and the horseshoe 
orbit. In the middle graph the typical behaviour for a body in a 
horseshoe orbit is shown which we can describe at the best in 
a rotating coordinate system (comp. Eig.[T]): close to the turn- 
ing point of the orbit the Trojan finds itself in the vicinity of the 
Earth, the semimajor axes is larger than the one of the Earth and 
then - when the libation angle changes from about 10° to 350° 
it reaches values smaller than aEarth \ now the Trojans is on the 
other side of the Earth. In Fig. [121 in the lowest graphs an un- 
stable orbit is shown which - after leaving to be in a horseshoe 
orbit - is chaotic with librations between 0° to 360°. 




mrnm - 



Time in thousand years 
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Fig. 11. Libration amplitudes - defined as half of the total li- 
bration angle - in the region close to L4; the color stands for 
the largeness of the libration. 

To show the diff'erent kinds of orbits we depict three diff'er- 
ent ones for an inclination of the Trojan's orbit / = 1°: a tad- 
pole orbit deep inside the stable region (Fig. [T2l upper graphs), 
a horseshoe orbit close to the border of stability (Fig.[T2lmiddle 
graphs) and an escaping orbit outside the stable region after ini- 



Fig. 12. The temporal evolution of three orbits: a tadpole (up- 
per graphs) a horseshoe (middle graphs) and an escaping or- 
bit (lower graphs) In six panels we show the evolutions of 
the semimajor axis (left graphs) and the libration angles (right 
graphs) during 5000 years. 



4.3. Secular resonances 

In this part we present our preliminary results about the secular 
resonances to understand the structure of the stability diagram 
(Fig.©. A secular resonance happens when the precession rate 
of an object (m or Q) equals to one of the eigenfrequencies of 
the system (Murray & Dermott 1999) . In a secular resonance, 
the eccentricity or the inclination shows long-term oscillations. 
The eigenfrequencies of the system can be calculated through 
a linear analysis (Bretagnon 1974 ) or can be numerically deter- 
mined dLaskar I990I) . 

The reality is more complicated than the model adapted in 
the linear analysis. The secular perturbations among the plan- 
ets modify the precession rates. Due to their large masses the 
ones of the big planets are almost constant; on the contrary, 
the precession rates of the inner planets are strongly influenced 
particularily by Jupiter and Saturn. In Fig. [131 we show the dif- 
ferences between Jupiter's perihelion longitude (^175, as usually 
denoted) and the ones of Venus (^2) and the Earth (^173) in our 



Dvorak et al.: The Earth Trojan asteroid 2010 TK7 



9 



simulations. Obviously, during the periods ~ 0.05 - 0.35 Myr 
and ~ 0.45-0.75 Myr, tu^ and tu^ have nearly the same preces- 
sion rate as Jupiter. This makes the secular evolution of aster- 
oids in this planetary system complicated, as we will see below. 



360 
300 
240 
fe"" 180 
120 
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0.0 0.2 0.4 0.6 0.8 1.0 

Time (Myr) 

Fig. 13. The differences between perihelion longitudes of inner 
planets and of Jupiter. The abscissa is time and the ordinate is 
TJTk-TJTs, with ^ = 2 for Venus (thin black curve) and ^ = 3 for 
the Earth (thick red curve). 

At high inclinations the Kozai resonance causes large 
coupled variations of the eccentricity and the inclination 
(IBrasser et al. 20041) . which finally results in instability. At 
low and moderate inclinations the secular reso nances and 
y4 were found to afi 'ect the Earth Trojan^ dMorais 2OOII: 



Brasser & Letho 2002h . We notice there are some disagree- 
ments in the locations of these secular resonances between dif- 
ferent papers. To depict a detail resonance map we would like 
to leave it to a separate paper. Here we only illustrate the evo- 
lution of four typical orbits to show the eff'ects of secular res- 
onances that contribute to form the structures in the stability 
diagrams. In Fig. [8] these are the four orbits with semi-major 
axis a = 0.9995 AU and inclination / = 10°, 22°, 35°, 42°; the 
evolutions of the resonant angles are presented in Figs. [14] and 
[T5I respectively. The signal of the eccentricities can be found 
in Fig. [161 Obviously, the orbits with / = 10° and / = 35° are in 
the stable windows, while the orbits with / = 22° and / = 42° 
are in the unstable gaps in Fig. [8] 

Our preliminary calculations show that the most important 
secular resonances aff'ecting the Earth Trojans are y2, V3, y4 and 
vs . According to the linear theory of secular perturbation, the 
averaged temporal variation of the eccentricity of a celestial 
body perturbed by a perturbing planet can be roughly approxi- 
mated by 



< e >= C sin(A^ir^) = C sin(m -TUi) , 



(6) 



where C is a negative constant deter mined by the mass of the 
planet and the Laplace coefficients (iMurray & Dermott 19991: 




0.3 0.4 
Time (Myr) 

Fig. 14. The evolution of apsidal diff'erences between the 
Trojans and the planets: /^tu2 - tu-tu^ (black), I^tu^ - tu-tj^ 
(red), Ati74 - tu - tu/^ (blue) and kru^ - tu - tu^ (green). The 
upper panel is for / = 10° and the lower panel is for / = 35°. 



^ A Trojan precesses at the rate of the eigenfrequency gs (or ^4), 
which is the precession rate of the Earth (or Mars), is in the vs (or y^) 
secular resonance. 



iLi et al. 2006h . tu and vii are the perihelion longitudes of the 
celestial body and the planet, and b^vji - tu - rui with / = 
1,...,5. If 0° < b^TUi < 180° then < ^ > is less zero and if 
180° < I^TUi < 360° then < ^ > is greater zero. 

For the orbit with / = 10°, from t = 0.05 Myr to 0.3 Myr the 
angles Aiir2, I^tu^ and I^tu^ librate, all around values ~ 40° (see 
Fig. [141 upper panel). Such libration aff'ects (but not excites, 
according to Eq. O the eccentricity of the Trojan. The over- 
lap of these secular resonances introduces chaos to the motion, 
which can be seen in the irregular eccentricity evolution (black 
in Fig. [161). After r = 0.3 Myr the and vs vanishes as the 
corresponding resonant angles circulate, but I^tua, begins to li- 
brate around a value larger than 180°, i.e. y4 resonance appears. 
Although this y4 resonance may increase the eccentricity since 
l^TU^ > 180°, this eccentricity pumping eff'ect is off'set by the vs 
resonance, which is continuously present with the Am3 < 180°. 

Similar evolution happens for / = 35° (see Fig. [141 lower 
panel). In this case, before ^ = 0.1 Myr however, vs and vs li- 
brate around high values, resulting in eccentricity excitation in 
this period (blue in Fig. [161). The y4 resonance does not happen 
at all, while vs and vs are absent for about 0.1 Myrs from 0.4 
to 0.5 Myr. Most of time, the resonant angles are smaller than 
1 80° and the eccentricity stays low. 
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Time (Myr) 

Fig. 15. The evolution of apsidal differences between the 
Trojans and the planets: ATjj2 (black), Atjjs (red), ATjj4 (blue) 
and AuFs (green). The upper panel is for / = 22° and the lower 
panel is for / = 42°. 



The situation is quite different for / = 22° (see Fig. [TSl 
upper panel). The y4 resonance exists with the resonant angle 
around 180° all the way before the Trojan is expelled from the 
1:1 MMR at ^ ~ 0.63 Myr. No other secular resonances are 
present until t ~ 0.35 Myr and the Trojan's eccentricity varies 
in correspondence with the circulation of the main secular an- 
gles ATJTk (k = 2, 3, 5). The critical eccentricity excitation hap- 
pens when the vs resonance sets in at ^ ~ 0.35 Myr. In this secu- 
lar resonance the eccentricity reaches 0.25 at ^ = 0.55 Myr (red 
in Fig. [16]), where the Trojan may approach to Venus closely. 
And this destabilizes the Trojan orbit. 

For / = 42° (see Fig. \T5\ lower panel) the resonant an- 
gles of y2,V3 and vs librate around high values, resulting in 
a quick excitation of the Trojan's eccentricity to 0.45 (red in 
Fig. [16]). Such a high eccentricity makes its orbit to cross the 
one of Venus, the one of Mars and finally the Trojan becomes 
unstable. 

Based on these typical orbital evolutions we argue that the 
secular resonances V2, V3, y4, vs are all involved in determining 
the Trojans' stability. Among them the vs and vs play more 
important roles than the others. 
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Fig. 16. The eccentricity evolution for four orbits. The solid 
curves represent the orbits with initial inclination / = 10° (solid 
black), 22° (dash red), 35° (dot blue) and 42° (dash-dot green), 
respectively. The vertical dash and dot lines mark the moments 
when the orbits escape from the 1:1 MMR for / - 22° and 
/ = 42°, respectively. 

5. Conclusions 

We investigated in full detail the stability of the recently found 
Earth Trojan asteroid 2010 TK7 in different dynamical models. 
We derived a symplectic mapping based on the spatial elliptic 
restricted three-body problem to see the location of the mean 
orbit in the context of the phase space structure. We tested the 
validity of the simplified model with numerical integrations and 
found good agreement of the model on short time scales. We 
included the influence of the additional planets as well as the 
Moon to see their influences on the mean orbit of 2010 TK7. 
Next, we performed a detailed numerically study to propagate 
the orbit of the Trojan Earth asteroid both, forward and back- 
wards in time together with 100 clone orbits to incorporate the 
possible errors of the orbital parameters induced by the obser- 
vations. From the detailed study we are able to state the prob- 
ability of capture and escape of the Earth Trojan together with 
the estimation of the time the asteroid will stay close to the 
Lagrangian point L4 of the Earth. To this end we investigated 
in full detail not only the regime of parameters and initial con- 
ditions close to the actual found asteroid but rather on a grid of 
initial conditions and in the parameter space (a, /) using thou- 
sands of simulations, based on fictitious Earth-Trojans. With 
this we were able to define the region of stability and instabil- 
ity on both short and long time scales. 

The main results of the present study can be summarized 
as follows: we can confirm the result of (Connors et al. 201 1) , 
the asteroid 2010 TK7 lies in the tadpole regime of the Sun- 
Earth system. The orbital parameters indicate that most prob- 
ably the asteroid became an L4 Trojan some 1,800 years ago 
and it will jump to the L5 neighbourhood or a horseshoe orbit 
in about 15,000 years in future. Before it moved to the cur- 
rent tadpole orbit this asteroid was captured into the 1:1 MMR 
around 30,000 years ago, and it may stay in the resonance for 
another 200,000 years. The total life-time of the asteroid (be- 
ing in the 1:1 MMR) is less than 0.25 Myr. As a Near-Earth- 
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Asteroid, the closest approach of this asteroid to the Earth will 
be larger than 70 times the Earth-Moon distance when it is on 
a Trojan-like orbit. On short time scales it is possible to predict 
the orbit in the spatial elliptic restricted three-body problem 
(Sun-Earth-Asteroid) on intermediate time scales the influence 
of the Moon has to be taken into account. In addition, on longer 
time scales the influence of the other planets of our Solar sys- 
tem cannot be neglected, since the asteroid is on a chaotic orbit 
with close encounters to the unstable equilibrium L3 . We ex- 
pect the discovery of further interesting objects in the vicinity 
of the L4 or L5 equilibria of the Sun-Earth system. The long 
term integrations show: for low inclinations the stability region 
extends in terms of semi-major axis up to +.01 AU. The size 
of it decreases with increasing inclination up to a threshold at 
about / = 20°. Another stability window opens up between 28° 
to 40° and disappears soon after. The preliminary results of our 
frequency analysis indicate that the y2, V3, y4 and vs secular res- 
onances are deeply involved in the motion of the Earth Trojans, 
and we show this by several examples. 

From a qualitative point of view our simulations indicate 
that 2010 TK7 is situated inside an unstable region where all 
sorts of orbits like horseshoes, tadpoles and jumping ones are 
possible. We show - this is well confirmed - the Earth Trojan 
moves on a temporary captured orbit. On the other hand it is 
unclear why there are no other Trojans of the Earth found up to 
now because there exist large stable regions for small inclina- 
tions, but also for larger inclinations Trojans may exist for very 
long time scales. 
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